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Abstract 

We report on the influence of submerged bubble clouds on the remote sensing properties 
of water. We show that the optical effect of bubbles on radiative transfer and on the estimate 
of the ocean color is significant. We present a global map of the volume fraction of air in 
water derived from daily wind speed data. This map, together with the parameterization of the 
microphysical properties, shows the possible significance of bubble clouds on the albedo of 
incoming solar energy. 
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1 Introduction 

"The effect of bubbles on the color of the sea may be observed in breaking waves... Where a 
great many bubbles have been entrained by a breaking wave it is white. But where there are 
fewer of them it is blue-green or green, brighter than the sea but not as bright as the foamiest 
parts of the wave. Even after a wave has broken and the water is again quiescent, a pastel 
green patch often remains, slowly fading into the surrounding sea as the bubbles dissipate. 
Thus the effect of bubbles on the color of the sea is similar to that of solid particles (B ohren, 



1987)." Bubbles within the water and foam on its surface (Bukata, 1995; Frouin et al., 1996; 
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Stramski, 1994) can predominate in determining the radiative transfer properties of the sea sur- 



face at higher wind speeds. However, there is a limited knowledge about the radiative transfer 
properties of bubble clouds, their inherent optical properties (IOP), and their global climatol- 
ogy. Mobley (1994 ) and Bukata (1995 ) discuss qualitatively the surface properties of bubble 
clouds. Frouin et al. (1996 ) performed spectral reflectance measurements of sea foam at the 
Scripps Institution of Oceanography Pier. They observed a decrease of the foam reflectance 
in the near-infrared and proposed that the foam reflectance can not be decoupled from the ref- 
electance by bubbles. [Stramski (1994 ) concentrates on light scattering by submerged bubbles 
in quiescent seas and shows the scattering coefficient and the backscattering coefficient at 550 
nm in comparison with scattering and backscattering coefficients of sea water as estimated 
from the chlorophyll-based bio-optical models for Case 1 waters. In this exploratory paper, we 
report on the influence of bubble clouds generated by breaking waves on the remote sensing re- 
flectance and calculate not only the inherent optical, but also apparent optical properties using 
the radiative transfer model. We show that the optical effects of bubbles on remote sensing of 
the ocean color are significant. Furthermore, we present a global map of volume fraction of air 
in water. This map, together with the parameterization of the microphysical properties, shows 
the significance of bubble clouds on the global albedo of incoming solar energy. By proxy, we 
show the influence of the bubble clouds on the remote sensing retrieval of organic and inor- 
ganic components of the natural waters. It is worth mentioning that the bubble clouds coincide 
with the upper range of the euphotic zone and will, therefore, contribute to the dynamics of 
the upper-ocean boundary layer, heat distribution, and sea surface temperature (Thorpe et al., 
1992). In fact, our initial motivation for this work was an observation that the asymptotic 
radiance distribution is established close to the ocean surface in apparent contradiction with 
theoretical studies ( flatau et al., 1999 ). Thus, the light field must become diffuse at shallower 
depths than usually modeled. This leads to search for alternative mechanisms influencing the 
light distribution. Thus, the importance to light scattering of the bubble clouds goes beyond 
the remote sensing issues considered in this work. 

In the next section, we discuss in more detail the microphysical and morphological prop- 
erties of bubble clouds, because they have a direct bearing on their optical properties and 
radiative transfer. 



2 Physical properties of bubble clouds 



2.1 Morphology of bubble clouds in natural waters 

Indi vidual bubble clouds are generated by breaking waves, persist for several minutes ( Thorpe, 
1995), and reach to mean depths of about AH S , where H s is the significant wave height, but 
with some clouds extending to about 6H S . There is evidence that at high wind speeds, sepa- 
rate bubble clouds near the surface coalesce, producing a stratus layer (Thorpe, 1995). Fig [j] 



is based on a sonograph of [Thorpe (1984| ). The "bubble-stratocumulus" (b-Sc) is often ob- 
served by acoustic means ( farmer and Lemon, 1984 ; [Thorpe, 1995 ). The depth of the b-Sc 
layer is related to the wind speed and wind variability, but more specifically it is set by larger 
waves, such us those breaking predominantly in groups ( Thorpe, 1995| ). The "stratus layer" 
description should not be taken too literally. For sufficiently high winds there will be signifi- 
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cant concentrations through out the upper layer, but the variability within this layer can be very 
high. 

Wind Velocity B m/s 




Figure 1: Bubble-stratocumulus (b-Sc) layer generated by breaking waves with inhomogeneous 
deeper intrusions. The depth of the cloud is related to the maximum significant wave height. 



2.2 Bubble cloud climatology 



According to [Thorpe et al. (1992| ) in the absence of precipitation and in wind speeds exceeding 
about 3 m s _1 , wave breaking generally provides the dominant source of bubbles. The wind 
speed Wio at 10 m above the mean sea surface level is used to parameterize the volume frac- 
tion of air in water, / = V a i r /V w , where V a \ v is the volume of air, and V w is the volume of 
water. This parameterization is an approximation of a more complex wind-wave relationship. 
Figure || shows the volume fraction of air in water estimated from the W\o winds for January of 
1992. The volume fraction and the wind speed are assumed to follow a non-linear relationship 
flWalsh and Mulhearn, 1987| ) 

/ = /o< 4 + /i CD 

The coefficients /o and f\ were calculated at 2m assuming that f= 10~ 8 for W\q = 6.2m s _1 
and f= 1(T 7 for W w = 10.5m s" 1 flWalsh and Mulhearn, 1987|). |Vagle and Farmer (1992[) 



show that the volume fraction decreases with depth, changing from about 10 at 0.3 m to 
10~ 7 at 2.7 m. We base our parameterization on these findings and extrapolate ([[]) to near-sea- 
surface depth using exponential fit. 

The monthly averages of / were obtained by employing 1992 daily surface winds from 
the NCAR/NCEP ( Kalnay et al., 1996| ) reanalysis project, and averaging the daily volume 



fractions for each month. Thus, Figure ^| is based on the variability of the wind field on the 
scale of one day. In the winter of the northern hemisphere, one can observe maxima associated 
with the midlatitude storm tracks in the Northern Pacific. Cyclogenesis, common in western 
parts of the oceans during the winter, contributes to mixing and a large bubble cloud volume 
fraction. This can be observed to the east of the North American continent. The Intertropical 
Convergence Zone (ITCZ) region, with its associated deep convection, may also be a region 
of enhanced production of bubbles. The winds of the Southern Ocean have a strong effect on 
bubble formation during both summer and winter. In subtropical regions, to the west of the 
continents, the subsidence associated with the descending branch of the Hadley circulation is 
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responsible for the relative minimum in /. It should be stressed that these results are qualitative 
and that they c an be improved by more detailed breaking wave climatology models ( Kraus and 
Businger, 1994). The data in Fig |2| is indicative of regions where bubble clouds are potentially 
important in the interpretation of remotely sensed reflectance. 



2.3 Optical thickness 

The size distribution of bubble clouds determines the optical thickness and is, therefore, one of 
the most critical parameters entering the theory. There are two assumptions which simplify the 
development here: (a) we consider light scattering in the geometric optics regime for which the 
ratio proportional to bubble radius to wavelength x = 2irr / A is large and (b) we assume that 
there is an effective radius r e g which determines the optical properties of the size distribution. 
Both (a) and (b) are quite probable. The size parameter, x = 50, corresponding to a bubble 
radius of approximately 5 fim, is already in the geometric optics regime, and r cr r = 10/mi 
will satisfy (a). For bubbles with an effective radius of r e g, the volume attenuation (equal to 
scattering for a non-absorbing sphere) can be expressed as 

, n N Kir N of s 

b = QscaT^S = 2 — — s = 2/-. (2) 
Vw Vw Vair V 

Thus 

f 

6 = 2— (3) 

where s = irr 2 is the cross-section of a bubble with radius r, v is the volume of such a bubble. 
N is the number of bubbles in volume V w of water and / is the fraction of air in a volume of 



water. The effective radius is defined as r c g = v/s ( |Stephens et al., 1990| ; [King et al., 1 993| ; 



Bricaud and Morel, 1986D - Scattering efficiency defines how much of incoming light is being 



"blocked" by a particle by scattering processes, for large size parameters Q sca , tends to 2. This 



issue is discussed in detail by Bohren and Huffman (1983) 



The physical significance of / comes from the fact that it is determined by the large scale 
forcing such as the wind field. Thus, for a given synoptic or climatological setting, the mixing 
ratio / is, to some extent, pre-determined. On the other hand, the effective radius r e g depends 
on processes of much smaller scale then the large scale. These processes are coagulation, 
coalescence, coating by organic material, saturation, buoyancy, pressure, etc. Thus, Eq. ^de- 
fines the optical properties of bubble clouds on both the large- and sub-scales. The expression 
b = 2/ /r e g holds for poly dispersions, and the only difference with a monodispersion is that 
r e ff is defined via distribution averaged s and v. 



2.4 Effective radius of the size distribution 

Numerous observations of bubble size distributions are reported in the literature based on 
acoustic, photographic, optical, and holographic methods ( Wu, 1988^). A kulichev and Bu- 
lanov (1987) summarize results from 22 experiments using different techniques. As the ori- 
gins of bubbles are biological (within the volume and at the bottom), as well as physical (at the 
rough surface), we may expect large regional and temporal variations of bubble concentration 
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between coastal and open oceanic waters and between plankton bloom or no-bloom conditions 
( [Thorpe et al., 1992Q . 

Figure || presents a comparison of bubble spectra under breaking waves and quiescent sea. 
Recently published observations, using laser holography near the ocean surface, have shown 
that the densities of 10 to 15/im radius bubbles can be as high as 10 6 (per cubic meter per 
micron radius increment) within 3 m of the surface of quiescent seas ( O'Hern et al., 1988 ). 
These results are plotted as solid squares connected with a vertical solid line. The majority of 
bubbles injected into the surface layers of natural waters is unstable, either dissolving due to 
enhanced surface tension and hydrostatic pressures or rising to the air-water interface where 
the bubbles break ( Johnson and Wangersky, 1987j ). However, bubbles with long residence 
times, i.e. stable microbubbles, have been observed. For example Medwin (1977 ) observed 
nearly 2.5 x 10 6 of bubbles per cubic meter in the radius range 18 — 355 /im for small wind 
speeds. One of the stabilization mechanisms ( Mulhearn, 1981 ; Johnson and Wangersky, 1987 ) 
assumes that the surfactant material is a natural degradation product of chlorophyll, present 
in almost all photosynthesizing algae. Isao et al. (1990 ) have observed very large populations 
of neutrally buoyant particles with radii between 0.1 — lfim. Johnson and Wangersky (1987 ) 
and [Thorpe et al. (1992| ) proposed another stabilization mechanism based on monolayers of 
adsorbed particles. Numerical modeling ( Thorpe et al., 1992 ) can be used to study the effects 
of water temperature, dissolved gas saturation levels, and particulate concentrations on the size 
distribution of subsurface bubbles. The results of such numerical models provide additional 
evidence for the existence of a small size bubble fraction which is not adequately measured by 
acoustic or photographic techniques. The dashed line on Fig. || presents the mean concentra- 
tion in the model steady-state for the water temperature of 0C. It can be seen that maximum 
concentration is around 15/um and it is 2.5 orders of magnitude larger than for lOO^m bubbles. 
Other results presented (Figure |3j> are those of Johnson and Cooke (1979 ) observations at 4m 
in wind speed ll-13m s _1 (open squares). In situ acoustic measurements of microbubbles 
at sea by pVledwin (1977 ) are plotted as solid triangles. The solid triangles joined by a solid 
line are concentrations at 4m depth and 3.3ms -1 wind speed. These spectra were obtained on 
August 7, 1975 in Monterey Bay. The solid triangles are for midafternoon, February 10-16, 
1965 at Mission Bay, San Diego, 3m below the surface, and in 1.7-2.8m s~ x winds. The open 
triangles are from Baldy (1988 ) and include data at 30cm depth with wind and swell and at 
25cm with wind only. They are based on extensive laboratory experiments. The solid hexagons 
are data from Medwin and Breitz (1989 ) acoustic measurements obtained in the open sea at 
25cm depth under the water surface during 12m s _1 winds under spilling breakers. The lit- 
erature reviewed here and encapsulated in Figure || is a mixture of descriptions of the effect 
of active wave breaking, and of stabilised microbubbles observed largely in coastal situations. 
Currently, it is not clear how to parameterise the stabilized bubbles. From results such as those 
presented in Fig. [| in the case of "transient," open ocean bubbles, it can be estimated that 
the size distribution follows the power law dependence n(r) oc r~ a and a 4 (Walsh and 
Mulhearn, 1987; |Wu, 1988^ )" Even though small microbubbles may not contribute to the total 
mass, they may be important for the light scattering. Therefore, it is of interest to estimate 
the contribution of small bubbles to the optical thickness. Assuming that microbubbles are 
spherical (V^ = 4/37rr 3 iV) we can show that the optical thickness (r = bh) of a layer with 
geometrical thickness h is 

b oc hV^N 1 / 3 (4) 
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Figure 3: Comparison of bubble spectra under breaking waves and quiescent sea. The solid 
squares connected with the vertical solid line are from laser holographic data P'Hern et al. (1988"|). 



Open triangles are from [Baldy (1988D and include data at 30cm depth with wind and swell, and at 
25cm with wind only. The solid hexagons are data based on [Medwin and Breitz ( 1 989| ) obtained in 
the open sea at 25cm depth under the water surface during winds of 12m s _1 . The solid triangles 
joined by the solid line and solid triangles are average bubble densities measured under comparable 
conditions at different seasons. The solid triangles joined by solid lines are concentrations at 4m 
depth, 3.3m s _1 wind speed, obtained on August 7, 1975 in Monterey Bay. The solid triangles 
are for midafternoon, February 10-16, 1965 at Mission Bay, San Diego, 3m below the surface, 



1.7-2. 8m s— 1 winds. The dashed line is from a numerical model ( JThorpe et al., 19921J . The mean 
concentration in the steady-state model are plotted at temperature 0C. The open squares show 
( Johnson and Cooke, 1979Q observations at 4m in wind speed 1 l-13m s _1 . 



Both h and V a \ x are assumed to be fixed. It can be seen from Eq. ^| that the contribution to 
optical thickness by very small particles will be the same as that by very large particles if their 
concentration varies as N 3 (or steeper). Indeed is the case flWalsh and Mulhearn, 1987|; W u, 



1988a). 

What remains to be defined is the effective radius r c g. On the basis of measurements, 
estimates of small particle fraction, existence of background 1 /im microbubbles, modeling 
predictions, and steep slope of microbubble size distributions, we decided to use r c g = 10/rni 
as a typical "radiative response radius." This choice does not exclude existence of larger or 
smaller particles. The real value may be between l/xm for stabilized particles and 50/im for 
the open ocean and will depend on many environmental factors such as storm passage, wind 
speed, swell, wind variability, phytoplankton concentration, water temperature, gas saturation, 
and other properties. The 10-15 fold increase in the size of effective radius or similar decrease 
in the air volume fraction will reduce the importance of air bubbles to a very small effect. 
It may be instructive to calculate the scattering coefficient b for typical size distribution of 
bubbles in water. In such case we have 

b= / Qsc^r 2 — (5) 

or b = 2fs/v where s = J irr 2 dN(r) and v = / 4/37rr 3 dN(r), and dN(r) is the number of 
bubbles between r and r + dr in a volume of water V w . For typical size distribution of bubbles 
in water dN(r)/dr oc 1/r 4 we have 

r eS = Un(n/ro)/(---) (6) 

Consider r\ = 150 and r = 10 micrometers. This gives r e g- ~ 3.6r which shows that the 
choice of small bubble cut-off is important for the bubbles' optical properties. However, the 
choice of this cutoff is non-trivial because the spectrum of small bubbles is not understood well 
at present. 

We close this section with some general comment about the effective radius. It is not 
a directly measurable quantity and, in essence, it defines how dispersed the given amount of 
mass is. Scattering of incoming solar radiation is sensitive to total projected surface rather than 
to total mass. For this reason the effective radius is commonly used in radiation calculations. 
However, it should be stressed that the effective radius is a semi-inherent optical property 
because it carries information not only about the size itself, but also about the orientation of 
particles, their morphology, coating, size distribution, or departure from a spherical shape. In 
addition, estimates of effective radius, as used in satellite remote sensing, often contain bias 
due to unrealistic assumptions about other optical properties such as optical thickness, leakage 
of photons due to horizontal transfer, wavelengths, or technique employed in retrieval. In that 
sense the effective radius is also used (or abused) as a semi-apparent optical property. 



3 Results 
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3.1 Numerical model 



The numerical radiative transfer model used in this study is a slightly modified version of the 
Hydrolight 3.0 code ( Mobley, 1994 ; Mobley et al., 1994| ). In brief, this model computes from 



first principles the radiance distribution within, and leaving, any plane-parallel water body. 
Input to the model consists of the absorbing and scattering properties of the water body, the 
nature of the wind-blown sea surface and of the bottom of the water column, and the sun and 
sky radiance incident on the sea surface. 

Pure sea water absorption and scattering coefficients are determined from the data of P ope 
and Fry (1997). 35 model wavebands were specified to cover the 400-700 nm region with a 
typical resolution of lOnm. The water column was specified as infinitely deep. Up to 62 depth 
layers, extending to 50 meters, were specified with a resolution of 0.5 m close to the surface. 
A clear sky was assumed, but the diffuse sky radiance was included. Three or four component 
systems were considered, consisting of pure water, particulates with or without bubble clouds, 
and dissolved organic matter. The spectral absorption of dissolved organic matter was defined 
as 

o(A) = a(A ) exp[-0.014(A - A )] (7) 

where a(Ao) = 0.1 m _1 , Ao = 440 nm ( |Bricaud et al., 198l"l ). The phase function of phyto- 
plankton was defined as an average of Petzold's clear ocean, coastal ocean, and turbid harbor 



cases (Mobley, 1994; Tyler, 1977). The bubble cloud phase function was calculated for the 



real relative refractive index m = 0.75 as an average for the size parameter range between 



x = 100 — 300 with a resolution Ax = 1 using Wiscombe (Fiedler-Ferrari et al., 1991 ) code. 
Figure || compares Petzold and bubble phase functions. The main reason for the size distri- 
bution average was to remove transient spikes. This normalized phase function is scaled by b 
calculated from Eq. |3[ The volume scattering and absorption coefficients for particulates were 
dete rmined from the Gordon-Morel model (^lobley, 1994 ; pordon and Morel, 1983|; M orel, 



1988) 



bp = ^o.3C - 62 (8) 
A 



a Uasei (A) = MA) + 0.06a*(A)C u - b5 ] [1 + 0.2 exp(-0.014(A -440))] (9) 

Here a w (X) is the absorption coefficient of pure water, a* is the chlorophyll-specific absorption 
coefficient, and C is the chlorophyll concentration in mg m -3 (Mobley, 1994). The chloro- 



phyll concentration was set as constant (well-mixed) with depth and equal to 0.8mg m~ 3 or 
0.08mg m~ 3 . We used Case 1 water parameterization but assumed that not all dissolved or- 
ganic material is correlated with the chlorophyll concentration. In an apparent contradiction, 
the wind speed which defined surface reflectance and transmittance functions due to the wind- 
blown water surface was set to 0. The reason for this was to estimate the effect on remote 
sensing properties of the sub-surface bubble clouds. However, we investigated the sensitivity 
of the reflectance to change in the wind speed between and 10 ms _1 , and the effect was 
small compared to the influence of bubbles. The azimuth direction was divided into 24 equally 
spaced sectors, the zenith-nadir range was divided into 20 equally spaced sectors. The profile 
of the bubble cloud volume fraction was determined by the expression 

f(z) = Hh/hf do) 
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Figure 4: Unpolarized phase function for a uniform distribution of bubbles between size parameter 
100 and 300 and refractive index m = 3/4 as a function of scattering angle 9. Also shown is 
Petzold's phase function. The bubble phase function has peak at 70 degrees. 
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where f = f(z ) = 10 6 , fx = f{zi) = 10 7 ,z = 0m,zi=8 m, and p = (z - z )/ [z\ - 
zo). Our choice of the speed of attenuation is perhaps unreasonably gradual, except for the 
case of stable microbubbles. On the other hand this choice is of secondary importance for the 
remote sensing properties of bubble clouds which are dominated by the surface volume frac- 



tion. In-water asymptotic light field is discussed by Flatau et al. (1999 ). The scattering was 
conservative. The volume attenuation was determined from the asymptotic expression (||). The 
probability distribution function for scattering was discretized and was set as constant through- 
out the depth range. We calculated the asymmetry parameter to be g m 0.85. This parameter 
defines the probability of photon scattering towards the forward or backward hemisphere. It is 
equal to 1 if all photons are scattered forward and to - 1 if all photons are scattered backward. 

3.2 Remote sensing reflectance 

The reference model runs were performed with the 3 component system of phytoplankton, 
pure water, and DOM with a constant chlorophyll profile. Two cases were computed (but only 
one is presented), representative of clear coastal (C = 0.8) and oceanic (C = 0.08) water. The 
remote sensed reflectance is defined as 

where is the downwelling irradiance onto the sea surface, and L w {\) is the upwelling 



water-leaving radiance. The remote-sensing radiance is a measure (Mobley, 1994) of how 
much of the downwelling light that is incident onto the water surface is returned into the zenith 
direction. The remote sensing reflectances are plotted in Fig. |[ The physical importance of 
the remote-sensing reflectance is evident from asymptotic theories which relate i? rs to inherent 



optical properties (Zaneveld, 1995) 



P(n - e m ) 

K TS oc (12) 

a 

where is the phase function, 6 m is related to the sun zenith angle, and a is the absorp- 
tion coefficient. Thus, R rs is approximately proportional to the probability of back- or side- 
scattering, and inversly proportional to the absorption of water column. Figure || shows the 
remote-sensing reflectance for the 3- and 4-component system with and without ocean bubbles 
but with constant pigment amount. The total single scattering albedo is strongly influenced 
by scattering from air bubbles. This leads to enhanced reflectance at all wavelengths. The 
results for C = 0.08 (not presented) show even larger sensitivity. In Fig || the gray rectangles 
indicate bands (wavelenghts) which are used by the current ocean color satellite instrument 
(SeaWiFS). It is of interest to comment on the remote sensing of pigments and bubble cloud 
retrievals. Consider an algorithm based on the ratio of remote-sensing reflectance and define 
the ratio of remote-sensing reflectances without (Chi) and with (Chl+b) microbubbles as 

ratio(A) = Rg hl /Rg hl+b . (13) 

Figure || shows ratio(A). Performance of the pigment algorithms based on the ratio of re- 
flectances will depend on ratio(Ai)/ratio(A2). Submerged microbubble clouds seem to be 
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wavelength-selective and even the ratio algorithms may require slight systematic correction. 
Given the increased sensitivity of the current generation of ocean color instruments, the ab- 
solute value of the radiances at the top of the atmosphere can be used for pigment retrievals. 



4 Summary 

Our calculations indicate that the optical effects of submerged microbubbles on the remote 
sensing reflectance of the ocean color are significant. These results are of importance for the 
retrievals of pigments from the ocean color measurements and for studies of the energetics of 
the ocean mixed layer. We provide information on how to reduce the systematic error due to 
microbubbles in pigment retrieval schemes via the ratio (A). We also derive apparent optical 
property of bubbles - remote sensing reflectance - for the whole solar spectrum. This AOP 
is directly observable by the satellites and remote sensors. We expect that these and similar 
AOPs will have to be invoked in the case of hyperspectral retrievals for Case 2 waters where the 
signals from minerals, bubbles, chlorophyll, and dissolved organic material (CDOM) are not 
well correlated. New algorithms for current satellite instruments such as MODIS and SeaWIFS 
should employ this information. 

We also present a global map of the volume fraction of air in water derived from daily wind 
speed data. We expect that such a map can be improved by knowing the day-to-day variability 
of the wind-wave relationship and better estimates of the volume fraction. 

The paper is exploratory. Therefore, it is perhaps worth playing advocatus diaboli and 
speculate why the bubble clouds may not be important, at least in current satellite ocean color 
retrieval practice. Here are some reasons: (1) The high wind and clouds are correlated. This 
masks (bias) the effect of bubbles as observed from the satellites; (2) Both whitecaps (Wu, 
1988b; Gordon and Wang, 1994| ; prouin et al., 1996D and bubble clouds are correlated via their 
dependence on wind speed. Therefore, our results, as well as the hypothesis of Frouin et al. 
(1996), indicate that reflectance of foam has to be considered together with the reflectance due 
to bubble clouds. On the other hand, there are cases in which strong wind is not correlated with 
clouds. For example, the cross equatorial flow during the summer monsoon in the southern 
Indian Ocean is strong but the ITCZ position is in the northern hemisphere. 

It is interesting to note that stabilized, coated microbubbles are hypothesized to be corre- 
lated to phytoplankton and CDOM concentrations; we need parameterization of this process. 
The optical properties of the first several meters below the surface are difficult to measure 
and are often removed from data due to experimental problems such as ship shadow or wave 
activity. This is the region where more detailed studies are needed. 
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